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ABSTRACT 

We have analyzed Hubble Space Telescope spectroscopy of 24 nearby Active Galactic Nuclei (AGNs) to 
investigate spatially-resolved gas kinematics in the Narrow Line Region (NLR). These observations effectively 
isolate the nuclear line profiles on less than 100 parsec scales and are used to investigate the origin of the sub- 
stantial scatter between the widths of strong NLR lines and the stellar velocity dispersion cr* of the host galaxy, 
a quantity which relates with substantially less scatter to the mass of the central, supermassive black hole, and 
more generally characterize variations in the NLR velocity field with radius. We find that line widths measured 
with STIS at a range of spatial scales systematically underestimate both cr* and the line width measured from 
ground-based observations, although they do have comparably large scatter to the relation between ground- 
based NLR line width and cr*. There are no obvious trends in the residuals when compared with a range of host 
galaxy and nuclear properties. The widths and asymmetries of [O III] A5007 and [S II] AA6716,6731 as a func- 
tion of radius exhibit a wide range of behavior. Some of the most common phenomena are substantial width 
increases from the STIS to the large-scale, ground-based aperture and almost no change in line profile between 
the unresolved nuclear spectrum and ground-based measurements. We identify asymmetries in a surprisingly 
large fraction of low-ionization [S II] line profiles and several examples of substantial red asymmetries in both 
[O III] and [S II]. These results underscore the complexity of the circumnuclear material that constitutes the 
NLR and suggests that the scatter in the NLR width and cr* correlation can not be substantially reduced with a 
simple set of empirical relations. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
Seyfert 



1. INTRODUCTION 

The width of emission lines from the Narrow Line Re- 
gion (NLR) in Active Galactic Nuclei (AGNs) has recently 
received a great deal of interest because it may provide a 
reasonably accurate, although not precise, estimate of the 
host galaxy spheroid velocity dispersion cr*. The NLR a g = 
FWHM/2.354 may therefore be a reasonable 'tertiary' black 
hole mass M, estimator based on a series of empirical re- 
lations that originate with the M, - cr* relationship between 
the mass of a galaxy's central, supermassive black hole and 
the stellar velocity dispersion of the host galaxy 's spheroid 
JFerrarese & Merriti2000l Idebhardt et alJl2000ah and the re- 
lation between cr* and NLR line width from Nelson & Whittle 
( 1996, and earlier work by Whittle 1992b,c). Subsequent 
measurements and analysis have found that the M, - a* re- 
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lationship has an intrinsic scatter of no more than 0.3 dex in 
logM. and th at this up per limit is due to present me asurement 
uncertainties (Tremaine et al. 2002). Nelson (2000) used this 
relation and M, estimates from reverberation-mapping exper- 
iments to propose the use of a g measured fr om the \Q IUl 
A5007 line to estimate the black hole mass. iGreene & Hoi 
(2005) recently completed a direct comparison of cr* and 
the gas velocity dispersion a g for approximately 2000 AGNs 
from the Sloan Digital Sky Survey (SDSS) and showed that 
these widths are well correlated, although with considerable 
scatter. 

The empirical relation between a g and M, is important for 
a broad range of applications. These applications include es- 
timates of M. for local AGNs, for AGNs at high redshift, 
and the cosmic evolution of supermassive black holes. For 
luminous AGNs with substantial continuum emission, or for 
AGNs at high redshift, the widths of the narrow emission lines 
may be the only method to determine cr*. In conjunction with 
measurements of bolometric luminosity Lb„i, estimates of M, 
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can also be used to calculate the accretion luminosity of AGN 
in terms of the Eddington ratio Lboi/LEdd- At low-redshifts, 
black hole estimates from reverberation mapping have been 
shown to agree quite well with the slope of the slope of the 
M.-g, relation jGebhardt et alJ2000bHFerrarese et al.l2001l 
iNelson et alJl2Q04t lOnken et al.ll2004l) and have been used to 
calibrate an additional secondary (virial) M, estimator based 
on t he line w idth and luminosi ty of the b road, permitted lines 
(e.g. lKaspi et al]l2000ft . iBorosonl l l2003l) used this virial rela- 
tion to show that the [O III] FWHM could be used to estimate 
M, to within a factor of five. The virial and a g estimates of 
M, have been used to explore the black hole masses and ac- 
cretion rates in particular classes of AGNs. most notably the 
Narrow Line Seyfert Is (Grupe & Mathui 2004), as well as 
to sea rch for evolution in the M, - c* relation ( Shiel ds et alJ 
l200l . Measurements of evolution in the M, - a* relation, or 
skewness in the distribution of the scatter about this relation, 
could also provide valuable const raints on the for mation his- 
tory of supermassive black holes (Robertson et a ll2005t) . 

The utility of a g measurements as a proxy for M, in AGNs 
is largely limited by the substantial scatter in the relation be- 
tween er» and a g . The brightest line for a g measurements 
is the [O III] A5007 line, although this line does suffer from 
substantial blue-side asymmetries (e.g. Heckman et al. 1981) 
that dramaticall y af f ect th e quality of the a g - <r* correla- 
tion. Idreene & Hoi J2005I) found that a g /a* = 1.34 ±0.66 
for the full profile of the [O III] A5007, while after removal 
of the blue-side asymmetry the result for the line core is 
a g /a* = 1.00±0.35. The other bright NLR lines, [Oil] A3727 
and [S II] AA6716,6731, have lower ionization potentials and 
do not suffer from such substantial asymmetries. The scatter 
and quality of the <r g -<j* relation for these lines is comparable 
to the core of the [O III] li ne. While all o f these correlations do 
strongly confirm the Nelson & Whittle ( 1996) result that the 
kinematics of the NLR are dominated by gravity, the origin of 
the scatter is not clear. One potential origin of this scatter is a 
nonvirial contribution to the emission line widths. These con- 
tributions are clearly present in the blue-asymmetric [O III] 
lines, but because these asymmetries are significantly stronger 
in high ionization lines like [O III], compar ed to lower ion - 
ization lines like [S II] (e.g. de Ro bertis & Osterbrockll984l) . 
there is good evidence that a range of spatial scales contribute 
to the observed NLR kinematics. Beyond the gas kinematics 
that contribute to the [O III] blue wings, any nonvirial contri- 
bution to deviations from the a g - cr* correlation is evidently 
traced equally by the low ionization [S II] lines and the core 
of the higher ionization [O III] line, since these features show 
similar scatter. 

Other potential origins for the scatter in the a g - er* relation 
include the impact of compact radi o jets and tidal distortions 
of the host galaxies. Nels on & Whittle! fl 996) found that both 
of these quantities correlate with systematically larger a g at a 
given <T», although large-scale radio jets or tidal distortions are 
not present in most galaxies and therefore can not be responsi- 
ble for the bulk of the scatter. Small- scale radio jets ma y play 
a more significant role, however, as Ho & Peng (2001) find 
approximately 60% of Seyfert 1 galaxies qualify as radio loud 
when the nuclear radio and visible-w avelength flux is isolated 
from host galaxy emission. Whittle (1992c) showed that the 
[O III] lines are significantly broader in Seyferts with linear 
radio morphology and high radio luminosity, which suggests 
radio jets can affect the NLR velocity field. The spatial dis- 
tribution of the NLR gas within the bulge may also have a 
substantial impact on the observed value of a g . While the 



NLR gas must approximately trace the spheroid kinematics 
to produce this correlation, spatially-resolved images indicate 
that the N LR is oft e n approx imately confined to a plane (e.g. 
Pogge et al.ll989bl) . IWhittlel(ll992bl) found that rotation does 
contribute to the width of the [O III ] emiss i on lin e, although it 
does not dominate. Finally, Greene & Ho (2005) investigated 
if host galaxy morphology, local environment, star formation 
rate, AGN luminosity, and Eddington ratio correlated with the 
observed scatter and only found evidence for systematically 
larger a g at higher Eddington ratio. While all of these investi- 
gations have determined that tidal disturbances, the presence 
of a radio source, rotation, and Eddington ratio contribute to 
the observed scatter, none of them dominate. It would be ex- 
tremely valuable if the origin of the scatter in the a g - c* cor- 
relation could be identified and substantially reduced, as then 
cr* (and M.) could be estimated more precisely for AGNs. 

Further information on the origin of the scatter in the a g - cr* 
correlation, as well as a more general understanding of the 
NLR kinematics, may be gained from a high spatial resolu- 
tion study. Nearly all studies of NLR kinematics to date have 
employed ground-based, single-aperture measurements with 
aperture sizes on the order of a few arcseconds. Such aper- 
ture sizes correspond to several hundred parsecs in projection 
and can include most of the NLR. There is good evidence, 
however, for the presence of stratification in the NLR, such as 
observations that higher -ionization and highe r critical den sity 
lines tend to be broader dde Robertis & O sterbroc kTl986l) . the 
presence of a blue, asymmetric wing on the higher-ionization 
[O III] line and not on other, lower-ionization lines such as 
[S II], and evidence for mode st increases in electron d ensity 
from recent HST observations ( Barthe^alj2^0J^ jShields_et alJ 
120051) . The properties of [O III! led lHeckman et al.l (119811) to 
propose that asymmetries in the NLR are due to rad ial outflow 
and wind models of the N LR (e.g. iKrolik & Vrtileklll984l 
Schiano 1986; Smith 1993) are supported by recent, spatially- 
resolved HST STIS k inematics of NGC 4151 (INelson et alJ 
2000; Das et al. 2005). These authors find evidence for both 
a rotational component and two, kinematically-distinct radial 
outflo w components that app ear to decelerate at larger scales. 
Crens haw & Kraemerld20 00) similarly studied the central 400 
parsecs of NGC 1068 and found evidence for radial outflow 
driven by either wind or radiation pressure in the nuclear re- 
gion, followed by deceleration on larger spatial scales. 

In the present paper we analyze spatially-resolved spec- 
troscopy of a large sample of nearby AGN from archival STIS 
observations and measure the NLR kinematics as a function 
of aperture size. We investigate the origin of the scatter in 
the a g - <7* correlation through a range of line width mea- 
surements at increasing aperture size, ranging from the lim- 
its of HST resolution to an aperture size typical of ground- 
based studies. We also study line asymmetries as a function 
of aperture size to investigate the spatial origin of line asym- 
metries. The sample properties and selection are described 
next in |j2]and the nonstandard aspects of the data processing 
in ^3] We describe the line width and asymmetry parameters 
in 50] along with the results of these measurements for dif- 
ferent spatial scales in ^5] The implications of this work for 
the a g -a* correlation and the kinematics of the NLR are then 
presented in ^6] We summarize our results in the final section. 

2. THE SAMPLE 

We have adopted the sample of INelson & Whittle! (1 1995! 
hereafter NW95) as the parent sample of this investigation of 
small-scale NLR kinematics because they have obtained high- 
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quality ground-based measurements of many NLR emission 
features, as well as stellar velocity dispersions, for a wide 
range of galaxy types. Their sampl e is largely draw n from 
the 140 Seyfert galaxies discussed in Whittle ( 1992a), where 
their primary selection criteria was for strong Mg b (~5 175A) 
and Ca II triplet (~8550A) absorption features and weak Fe II 
emission (which confuses Mg b). With moderate spectral res- 
olution (80-230 km s" 1 FWHM) they measured the redshift, 
stellar kinematics (V* and <t*), and a variety of characteristics 
for the [O III] A5007 and [S II] A9096 lines (see {gj. These 
high-quality ground-based data are a critical point of compar- 
ison for our study of the small-scale kinematics. 

We searched the HST archive and found 24 AGNs from the 
NW95 sample with moderate-resolution spectroscopy of ei- 
ther the [O III] or [S II] emission lines (G430M and G750M 
gratings, respectively) through either an 0.1" or 0.2" wide 
slit. The spectral resolution of the G750L and G430L grat- 
ings (~500 km s" 1 ) is not high enough to resolve line width 
changes in nearly all of these sources, while in contrast the in- 
strumental resolution is substantially higher for G430M and 
G750M. The 0.1" and 0.2" slit widths have a spectral res- 
olution of ~ 1.5 pixels FWHM for point sources and the 
G430M and G750M gratings have dispersions of 0.28 and 
0.56 A/pixel. The instrumental velocity resolution FWHM 
is therefore 25 and 40 km s" 1 , respectively. The resolution 
is approximately a factor of ^2.7 times larger for an unre- 
solved source that uniformly fills the 0.2" slit and a factor of 
~ 1.7 times larger for the 0.1" slit. The STIS plate scale is 
0.05078" per pixel and the slit length is 50". The details of 
our sample, including host galaxy properties, program ID, and 
observing mode, are provided in Tabled Of the 24 galaxies, 
only NGC 2110 has observations with both the G430M and 
G750M settings. 

3. DATA REDUCTION 

Our sample is drawn from many different STIS observing 
programs and thus the objects were obtained with many dif- 
ferent observing sequences. The reduction steps therefore 
vary from galaxy to galaxy depending upon the number of 
dithered exposures, the number of CR-split frames, and if the 
observations were obtained prior to the primary (Side 1) STIS 
electronics failure on 2001 May 16. We found that the re- 
duced spectra taken directly out of the CALSTIS pipeline are 
not of optimum quality for any of our galaxies. In this sec- 
tion we discuss our modifications to the CALSTIS pipeline 
for the various observing sequences employed to obtain these 
archival data. 

3.1. Multiple Dithered Exposures 

Twelve of our 24 galaxies (those obtained for propos- 
als 7361, 8055, and 9143) include multiple exposures of 
the nucleus and employed integer pixel dithering along the 
slit. These exposures had only one CR-split frame, and thus 
the CALSTIS pipeline reduction did not include cosmic ray 
correction. To account for the cosmic rays, we broke out 
of the CALSTIS pipeline after the dark-subtraction, bias- 
subtraction, and flat-fielding for each exposure (using the 
_flt.fits calibration file). After aligning the exposures, we 
used the Laplac ian Cosmic Ray i dentification IRAF routine 
(L.A. Cosmic; van Dokkum 2001) to create a bad pixel map 
for each 2-D spectrum. The pixels flagged as cosmic rays 
by L.A. Cosmic were then rejected when the multiple spectra 
were combined. L.A. Cosmic occasionally mistook spectral 



features for cosmic rays, leaving emission lines blank in the 
final combined image. In these cases we replaced the empty 
pixels with the corresponding values from the original _flt 
image. Once the cosmic-ray cleaning was complete, we re- 
inserted this combined, cosmic ray-rejected file into the CAL- 
STIS pipeline for the final wavelength and flux calibrations. 

3.2. Multiple CR-splits 

When multiple images were obtained at the same pointing, 
the exposures could be combined with the ocrre ject step 
in the CALSTIS pipeline to produce a cosmic ray corrected 
image. We have run the spectra from galaxies with multiple 
CR-splits listed in Table [Dup through this cosmic ray rejec- 
tion step in the CALSTIS pipeline. However, even with multi- 
ple CR-splits, the pipeline output most often is still peppered 
with cosmic rays. We thus used L.A. Cosmic to zap the re- 
maining rays in the ocrreject output (_crj.fits image) before 
re-inserting the spectra into the pipeline for the wavelength 
and flux calibrations. 

3.3. Side 2 STIS Electronics Correction 

The ten datasets for proposal 9143 were obtained after 2001 
May 16, when the Side 1 STIS electronics failed. The Side 2 
electronics do not have closed-loop temperature control of the 
CCD and this results in a dark rate that varies with tempera- 
ture. Thus the dark calibration image is a poor approxima- 
tion to the real dark rates for hot pixels and the outputs of the 
CALSTIS pipeline have many strongly negative pixels after 
the dark-subtraction step. Unfortunately, even the dithering 
scheme employed for this proposal did not easily remove the 
negative 'holes' created in the images. 

To correct for these negative pixels, we roughly followed 
the method outlined by Pogge et al. (2005, in prepara- 
tion). From t he co mbined, cosmic ray-rejected image for each 
galaxy (see 33. 1> . we created a list of pixels that fell more 
than 3-er below the background median, which we fed into a 
custom IRAF script that replaced each negative pixel with the 
mean of its surrounding 3 x 3-pixel box. We then re-inserted 
this corrected image into the CALSTIS pipeline. For images 
obtained prior to the Side 1 electronics failure, hot pixels that 
remained after the dark correction were eliminated as part of 
the cosmic ray cleaning step. 

4. DATA ANALYSIS 

Once these steps were complete, we extracted 1-D spec- 
tra of different aperture sizes from the 2-D spectra using the 
CALSTIS xld routine, which performs a geometric rectifi- 
cation and background subtraction. For each galaxy we ex- 
tracted spectra with apertures of integer pixel sizes from 2- 
50 centered on the nucleus, corresponding to angular sizes 
of 0.10" -2.54". We chose to integrate progressively wider 
apertures, rather than attempt to measure the differential line 
profiles as a function of radius, because one of the main goals 
of our study is to characterize the transition between NLR 
profiles on very small scales and the ground-based apertures 
used in previous work. 

Figure[2presents three 1-D spectra for each of the 24 galax- 
ies in the sample. The leftmost column of panels plots each 
galaxy's spectrum summed over an 0.2" aperture, which ap- 
proximately represents the nuclear emission unresolved by 
STIS. The middle column of panels plots the sum of a 1" 
aperture centered on the nucleus. The 1" aperture was chosen 
as a compromise between the resolution of previous ground- 
based studies and the angular extent over which emission lines 
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TABLE 1 

Observation and Sample Properties 



Name Class cz Dist. D Morph. \ogL XAGHz ID Mode CR-SPLIT Apertures Refs. Refs. Refs. 

(kms -1 ) (Mpc) (") (T) (W/Hz) (") Dist. D Morph. 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



MKN 270 


2 


3137 


38.1 


65.8 


-2 


21.51 


9143 


G750M 


1 


52X0.2 


1 


1 


3 


MKN573 


2 


5166 


68.9 


101.9 


-1 


21.56 


9143 


G750M 


1 


52X0.2 


4 


2 


3 


MKN 686 


2 


4251 


56.7 


90.8 


5 


22.51 


9143 


G750M 


1 


52X0.2 


4 


2 


1 


NGC 788 


2 


4061 


54.1 


77.3 





— 


9143 


G750M 


1 


52X0.2 


4 


2 


3 


NGC 1052 


2 


1307 


18.0 


77.3 


-5 


— 


7403 


G750M 


3 


52X0.2 


2 


1 


3 


NGC 1358 


2 


4008 


53.4 


119.7 





— 


9143 


G750M 


1 


52X0.2 


4 


1 


3 


NGC 1667 


2 


4547 


60.6 


53.5 


5 


— 


9143 


G750M 


1 


52X0.2 


2 


1 


3 


NGC 2110 


2 


2249 


30.0 


42.5 


-3 


— 


8610 


G750M 


6 


52X0.2 


4 


2 


3 


NGC 2273 


2 


1929 


27.6 


203.3 


5 


21.84 


9143 


G750M 


1 


52X0.2 


1 


1 


1 


NGC 303 1 


1.8 


-39 


1.4 


1439.3 


3 


20.81 


7351 


G750M 


2 


52X0.1 




1 


3 


NGC 3227 


1.5 


1024 


20.6 


119.7 


1 


21.84 


7403 


G750M 


2 


52X0.2 




1 


3 


NGC 3516 


1.5 


2503 


38.9 


137.5 


-2 


21.85 


8055 


G750M 


1 


52X0.2 




1 


3 


NGC 3982 


2 


1155 


17.0 


109.2 


3 


21.54 


7361 


G750M 


1 


52X0.2 




1 


3 


NGC 4051 


1.5 


622 


17.0 


361.5 


4 


21.35 


8228 


G750M 


1 


52X0.2 




1 


3 


NGC 4579 


1.9 


1334 


19.1 


378.6 


3 


21.34 


7403 


G750M 


3 


52X0.2 




1 


3 


NGC 5347 


2 


2295 


31.1 


101.9 


2 


21.07 


9143 


G750M 


1 


52X0.2 




1 


3 


NGC 5427 


2 


2733 


33.2 


181.2 


5 




9143 


G750M 


1 


52X0.2 




1 


3 


NGC 7682 


2 


5109 


68.1 


65.8 


2 


22.59 


9143 


G750M 


1 


52X0.2 


4 


1 


3 


MKN 348 


2 


4505 


60.1 


75.5 





22.17 


8253 


G430M 


2 


52X0.2 


4 


2 


3 


MKN 1066 


2 


3523 


47.0 


72.1 


-5 


23.16 


8253 


G430M 


2 


52X0.2 


4 


2 


3 


NGC 2110 


2 


2249 


30.0 


42.5 


-3 




8253 


G430M 


2 


52X0.2 


4 


2 


3 


NGC 4151 


1.5 


966 


20.3 


424.8 


2 


22.25 


8473 


G430M 


3 


52X0.2 


1 


1 


3 


NGC 5194 


2 


461 


7.7 


534.8 


4 


22.20 


7574 


G430M 


2 


52X0.2 


1 


1 


3 


NGC 5929 


2 


2250 


38.5 


65.8 


2 


22.22 


8253 


G430M 


2 


52X0.2 


1 


1 


3 


NGC 7674 


2 


8673 


115.6 


60.0 


4 


23.63 


8259 


G430M 


2 


52X0.2 


4 


1 


3 



REFERENCES. —Distances: (lt ^lTl^^: f2>|Wg^all ^TB; f3) liinlane& stall Eiffll : (4) H = 75 km s _1 Mpc -1 (adopted for consistency with Tidlv 19821). Size and Morphology: (1) Nilsou 1973i; (2) 
Voronlsov-Vclyaminov * Arkipov: 1974 i; (3)dc Vaucoidcurs el al 1991 1 

NOTE. — Col. (1): Galaxy name. Several common aliases are MKN 270 = N GC 5283. MKN 686 = NGC 5695, NGC 3031 = M81,and NGC 5194 = M51a. Col. (2): AGN type, where Type 1 AGN exhibit broad permitted lines and Type 
2 AGN have narrow permitted lines. Col. (3): [O III] profile median redshift frorr lNelson & WMttlelll99ll ([Offl] C80or [S III] C80 redshifts used for galaxies will] no [O III median redshift given). Col. (4): Distance. Col. (5): Angular 
diameter in arcsec of the 25 mag arcsee - ^ isophote. Col. (6): Morphological type (de Vaucouleurs numerical type). Col. (7): Radio luminosity (logjQ of the absolute spectral luminosity at 1.4GHz) from Condon ct al 2002 1. Col. (8): HST 
proposal ID for data used. Col. (9): STIS grating. Col. (10): Number of CR-splil exposures. Col. (11): STIS aperture. Cols. (12)-(14): References to published distances, sizes and morphologies. 



are detectable in most of the STIS spectra. Finally, the right- 
most column of panels shows the difference of the middle and 
leftmost panels. The difference spectra in these panels thus il- 
lustrate the line profiles outside of the nucleus and will be 
discussed below in the context of line widths and asymme- 
tries outside of the nuclear region. 

We applied a scalar throughput correction to the flux- 
calibrated 0.2" spectra (4-pixel extraction box height) pro- 
duced by xld because the default throughput correction table 
(referenced by the PCTAB keyword) does not include an entry 
for a 4-pixel extraction box height. When a 4-pixel extraction 
box height is requested, the throughput value for a 3-pixel box 
height is employed instead by default, although the through- 
put is approximately 10% larger for a 4-pixel box compared 
to a 3-pixel box. The default behavior of xld therefore over- 
estimates the flux in a 4-pixel box when it assumes the lower 
throughput appropriate for a 3-pixel box. To properly correct 
the 0.2" spectra, we interpolated the throughput values in the 
correction table to calculate the appropriate quantity for a 4- 
pixel box height for each instrument configuration and then 
applied this throughput correction to the xld output. Also, 
we note that while the throughput correction in the table was 
calculated from observations of a point source, the AGN in 
this sample are dominated by unresolved nuclear emission 
and this correction table is a reasonable approximation. This 
throughput correction is only relevant to the l"-0.2" and 0.2" 
spectra shown in Figure[2because the correction is negligible 
at 1". The correction also does not affect the width measure- 



ments because the spectra are only multipled by a constant. 

We have chosen to characterize the line profiles of the 
[O III] and [S II] emission lines as a function of aperture size 
through line widths measured at a range of heights, interper- 
centile velocity widths (IPVs), and a symmetry measures (e.g. 
iHeckman et al.ll981llWhittleil985al) . The definitions of these 
parameters, our measurement procedure, and their physical 
motivation, are discussed in detail below. The measurements 
of each of these quantities are listed in Tables 3 and 4 for 0.2" 
and 1" apertures. We present an illustration of these parame- 
ters in Figure |2 

4.1. AGN Continuum Subtraction 

The first step in measurement of the [O III] and [S II] line 
profiles is accurate subtraction of the continuum. This is par- 
ticularly important for characterization of any broad wings 
on the line profiles. In all cases, we defined continuum re- 
gions approximately 300-600 km s" 1 away from the lines un- 
der consideration. For all but four of our targets (NGC 3227, 
NGC 3516, NGC 4051, and NGC 4579), we found that the 
faint AGN continuum was well-approximated by a constant 
over this narrow wavelength range. For these four bright 
galaxies, we had sufficient signal-to-noise ratio in the con- 
tinuum to justify a power-law fit. The STIS slit size is suf- 
ficiently narrow that no significant stellar features appear to 
contaminate the continuum. 

4.2. Line Width Measurements 
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FIG. 1. — Profiles for the [S II] AA6717,6731 or [O III] A5007 lines for our sample are presented at aperture sizes of 0.2" (left column) and 1" (middle 
column), as well as for 1" —0.2" (right column). The best fit Gaussian profiles are shown as smooth, solid lines, while the individual components are shown 
as grey lines. Residuals are plotted below each fit. The vertical dashed line represents the position of the centroid of each line profile in the 1" aperture. A 
polynomial fit to the AGN continuum has been subtracted from NGC 3227, NGC 3516, NGC 4051, and NGC 4579. All other profiles have only had a constant 
subtracted. The profiles have been plotted in rest wavelengths using the cz values from in TablefTI Note that a Gaussian profile was not fit to NGC 7682 because 
the STIS G750M grating was centered too far to the blue and the broad, red wing of the [S II] A6731 line did not fall on the detector. 



We have characterized the mean velocity structure in the 
[O III] AA4959,5007 and [S II] AA6717,6731 emission lines 
through measurement of FW20, FWHM, and FW80, which 
are the line widths at 20%, 50%, and 80% of the peak height, 
respectively. These parameters, as illustrated in Figure [2J 



yield the base widt hs, core widths, and top profile widths 
described in Heck man et alJ ( 11981ft . We have linearly in- 
terpolated the flux between wavelength bins (^0.28A for 
the G750M and G430M gratings) to measure the widths. 
The value FWHM/2.354 was used for comparison with 
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FIG. 1. — Continued 



the stellar velocit y disp ersion measurements measured by 
lNelson& WhMel Il995). 

Several of the AGNs have sufficiently blended [S II] lines 
such that the flux between the 67 17 A, 6731 A doublet exceeds 
the 20 percent peak value for both lines. For the most severely 
blended of these, the FWHM value was undefined as well, 
and for NGC 1052 and NGC 2110 the FW80 value could 
not be measured. We find that the [S II] lines often become 
more blended with increasing aperture size and for several 



galaxies the FW20 and FWHM measurements become unde- 
fined above some aperture size. In these cases we measured 
the width values from a Gaussian fit to t he li ne profile (our 
Gaussian fitting scheme is described in 34. 3> . Tabled lists 
our width measurements of [S II] and/or [O III] for each of 
the sources at aperture sizes of 0.2" and 1.0" and indicates 
the values obtained from a Gaussian fit (with a g superscript). 
Figure|3]shows the value of these three quantities as a function 
of aperture size. 
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These width measurements are the observed widths and 
have not been corrected for the instrumental resolution. As 
noted in section |2] the FWHM velocity resolution of the 
G430M and G750M gratings are ~ 25 and 40 km s" 1 , respec- 
tively, for an unresolved point source and a factor of up to 
^2.5 times higher for a source that fills the slit. The most 
commonly used method to correct for broadening due to the 
instrumental resolution is to subtract the width of the instru- 
mental resolution in quadrature. This correction is can not 
be performed with great accuracy for these observations be- 
cause there is a substantial difference in resolution between 
unresolved and resolved emission and both sources of emis- 
sion likely contribute to the observed emission line profile. 
However, we can estimate the magnitude of this effect for the 
narrowest lines in our sample, which have FWHM ^120 km 
s -1 (see Table and represent the most affected measure- 



ments. If this emission were completely unresolved the cor- 
rection would be less than 5% and if it uniformly filled the 
slit it would be less than 30%. Because the two-dimensional 
spectra of these sources all indicate that the emission-line gas 
is centrally peaked, we estimate that the majority of the emis- 
sion in the smallest aperture is still unresolved. We therefore 
conclude that instrumental resolution makes at most a small 
contribution to our reported line widths. For the majority of 
our sample, the line widths are sufficiently broad that this cor- 
rection is less than 10% even with the most extreme assump- 
tion that the emission uniformly fills the aperture. We also 
note that the instrumental resolution correction is likely to be 
more important for the line c ore (FW80) tha n at smaller frac- 
tions of the peak height (e.g. Whit tlell9853) . 

In their study of 1749 AGNs spectra from the SDSS, 
iGreene & Hoi (|2005) found that the moments of Gaussian fits 
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FIG. 2.— A sample [O III] A 5007 profile (0.2" aperture) with definitions of the width parameters FW20, FWHM, FW80, IPV10, and IPV20, and the 
asymmetry parameters A10 and A20. Note the broad wings and blue asymmetry typically seen in AGN [O III] profiles. 



TABLE 2 
Line Width Measurements 



Nelson & Whittle ( 1 995 ) Aperture = 0.2" Aperture = 1 " 

Name ct* Line FWHM FW20 Line FW80 FWHM FW20 FW80 FWHM FW20 l"/0.2" NW95/1 

(km/s) (km/s) (km/s) (km/s) (km/s) (km/s) (km/s) (km/s) (km/s) (%) (%) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



MKN 270 


148 


[O III] 


396 


620 


[S II] 


100 


223 


457* 
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494* 


14 


56 


MKN 573 


123 


[S III] 


290 
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[S II] 


81 


176 
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91 
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420 


5.7 


56 


MKN 686 


144 


[O III] 


359 


609 


[S II] 


97 


192 


313 


99 


216 


370 


13 


66 


NGC 788 


140 


[O III] 


190 


295 


[S II] 


70 


140 


247 


77 


150 


246 


7.1 


27 


NGC 1052 


207 


[O III] 
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1085 


[S II] 
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966* 
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3.4 


26 


NGC 1358 
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410 
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361 
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NGC 1667 
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[0 III] 
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[S II] 
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462 
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[S II] 
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NOTE. — Col. (1): Galaxy name. Col. (2): Sldlar velocity dispersion measurement from NW95. Col. (3) Line used for width measurements in NW95. Cols. (4)-(5): FWHM and FW20 values from NW95. Col. (6): Line used for our width 
measurements. Cols. (7)-(9): Our FW80, FWHM, and FW20 measurements in an 0.2" aperture (the g superscript denotes measurements obtained from a Gaussian model of the line profile). Cols. (10M12): FW80, FWHM. and FW20 values 
measured in a \" aperture. Col. (13): The percent change in FWHM between the 0." 2 to \" aperture sizes. Col. (14): The percent change between our \ " FWHM measurement and theNW95 value. 
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FIG. 3. — Line width parameters FW20 (diamonds), FWHM (squares), and FW80 (crosses) are plotted against aperture size in units of physical distance 
(below) and angular size (above) for each galaxy in our sample. The NW95 measurements for FW20 (dashed), FWHM (dotted), and 2.354 X ct* (dash-dotted) 
are plotted as horizontal lines. Black symbols denote measured widths, while gray symbols denote values obtained from a Gaussian model of the line profile. 
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to the profiles provide more robust widths than direct mea- 
surements from the spectra. However, the SDSS data are typ- 
ically lower signal-to-noise ratio (SNR) and width measure- 
ments can be more uncertain in the low SNR regime. The 
STIS spectra used for our sample generally have sufficient 
SNR and resolution that we can confidently measure emission 
line widths without assuming any model for the line shape. 
Also, because the majority of the line profiles have a blue or 
redshifted wing and/or a broad central component, measur- 
ing widths directly from the data allows us to avoid approx- 
imating single values of the moments from a combination of 
multiple Gaussian components. 

The dominant uncertainty in direct measurement of the line 
widths is the continuum level because these lines hav e high 
signal-to-noise ratio. As noted above in Section |4~TI we fit 
the continuum by either a constant or a power law. The line 
width measurement uncertainty depends on both the uncer- 
tainty in the continuum level and the line profile shape, as for 
example a given continuum uncertainty will produce a smaller 
width uncertainty in a broad line than a narrow line. We esti- 
mate that the uncertainties in our direct width measurements 



are less than 5%. The width measurements that required 
Gaussian fits have a formal uncertainty of approximately 5% 
(Gree ne & Hol l2005). although this is only a true estimate if 
the lines can be correctly represented by Gaussians. Those 
cases where the fitting routine varied between direct width 
measurements and Gaussians (the switch between black and 
gray symbols in Figure [3} indicate that the uncertainties are 
at most 10%. The lack of substantial stellar continuum emis- 
sion in the extremely narrow STIS slit removes a potentially 
significant source of the uncertainty in width measurements 
based on ground-based observations, which may be particu- 
larly relevant for fainter emission lines such as [S II]. 

4.3. Area and Asymmetry Measurements 

Whittle ( 1985a) advocates the use of area measurements to 
define line width parameters, rather than simple cuts at vary- 
ing heights, because they have an integral nature and are thus 
smoothly defined and less sensitive to the presence of noise 
or the effects of instrument resolution. His definitions of the 
interpercentile velocity widths (IPV10 and IPV20) and asym- 
metries (A10 and A20) are also illustrated in Figure 2. The 
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median is the wavelength that denotes the center of area for 
the profile, and the lengths a, b, c, and d represent the separa- 
tion between the median and the profile's 10%, 90%, 20%, 
and 80% area values, respectively. The IPV10 parameter 
characterizes the base and win gs of the profile , much as did 
the FW20 width suggested by Heckm an et alJ (|T98 1 ), while 
the IPV20 parameter and higher percentage areas character- 
ize the line core. The A10 and A20 values serve to clearly 
quantify the profile's red or blue asymmetry. 

We have chosen to measure our area parameters from Gaus- 
sian fits to the l ine profi l es, rath er than directly from the data 
as discussed in Whittle (1985a). This choice was driven by 
the [S II] lines, which are mildly to severely blended in all 
of our sources. Because these lines are blended, we could 
no t measure the in terpercentile markers in the same manner 
as lWhittld dl985al) for the more isolated [O III] A5007 line. 
We experimented with taking the 10 and 20 percentile mark- 
ers of the total [S II] doublet area from the blue wing of the 
6717A line, and the 80 and 90 percentile markers from the 
red wing of the 6731 A line, and then scaling the widths a, 
b, c and d by the relative widths of the [S II] lines, but the 
potentially variable line ratio of the [S II] doublet made this 
approach uncertain. In the end, we found that the Gaussian 
profiles provided more robust measurements of the area pa- 
rameters, as we were able to isolate the separate contributions 
to the line profile from each of the [S II] lines. Also, the Gaus- 
sian fits allowed us to avoid the subtl eties in treating f ractional 
wavelength increments discussed bv lWhittld(IT985al) . 

Very few of these galaxies are well fit by a single Gaussian. 
The [O III] line often has an asymmetric blue wing in addi- 
tion to a broad central component, and so we allow for up to 
three Gaussian components in our fitting routines. With the 
exception of the sources with bright AGN continua, we fit a 
constant continuum term. We fit the [O III] AA4959 and 5007 
lines simultaneously with their wavelength separation and rel- 
ative strengths fixed to the theoretical ratio of 3:1 determined 
by atomic physics. For the [S II] doublet we also fit the two 
lines simultaneously and with a fixed separation, although we 
allow their relative line strengths to vary because the [S II] 
A6717/A6731 ratio is sensitive to electron density. In prin- 
ciple, the line widths of the 6717A and 6731 A lines may be 
different due to stratification in the NLR, so we initially al- 
lowed the widths of the two lines to vary in our fitting routine. 
However, we found that the widths of the two lines differed 
by less than 3%, which is within the error of the Gaussian fit 
parameters described in the previous section. We therefore 
fixed the widths of these two lines in the fitting routine. 

The fitting routin e employed here was adopted from 
iGreene & Hoi (|2005). The first step is to fit a single Gaus- 
sian to the profile, and then we allow a second central com- 
ponent centered between -5A and +5A of the line peak (al- 
though occasionally the component was best fit outside these 
limits). We keep this second component if the x 2 value for 
the fit improves by at least 20%. For [O III] we experimented 
with a third, blue component that was generally centered be- 
tween -20A and OA of the central component. Again we kept 
this component when x 2 improved by at least 20%. For the 
[S II] lines we only allowed for one additional component be- 
tween -10A and +10A of the line peak. The best fits from our 
Gaussian routine for each forbidden line profile are shown in 
Figure ^and the fit residuals are plotted below the spectrum 
in each panel. The IPV at 0."2 and 1" are provided in Tabled 
and plotted as a function of aperture size in Figure |4] while 



the asymmetry measurements are listed in Tableland shown 
in Figure |5] Unlike the width measurements, the IPV area- 
defined parame ters are re latively insensitive to instrumental 
resolution jWhitflell985al) . 

5. RESULTS 

All but three of the 24 galaxies in the sample have mea- 
surable emission outside of a central, 0.2" aperture, where 
0.2" corresponds to 8 - 115 parsecs for this sample (exclud- 
ing the very nearby NGC 3031). The three galaxies that lack 
detectable emission outside of 0.2" (MKN 686, NGC 3516, 
and NGC 5427) were observed with G750M and not G430M. 
The brighter [O III] line is always detectable outside of the 
central 0.2". The presence of significant emission outside of 
0.2" is illustrated in the rightmost columns of Figure^ which 
displays the difference between the sum of a 1" and 0.2" aper- 
ture centered on the nucleus. In addition to the three galaxies 
that do not exhibit detectable emission outside of 0.2", nine 
additional galaxies with [S II] observations have weak emis- 
sion on larger scales and six have substantial emission. How- 
ever, even emission lines that are quite weak outside of 0.2" 
can contribute to the line profiles in larger, integrated aper- 
tures. 

5.1. Line Widths vs. Aperture Size 

The line width and IPV values as a function of aperture size 
are shown in Figures [3] and |4] respectively. Figure [3] plots 
the radial dependence of FW20, FWHM, and FW80 for each 
galaxy, while Figure |4] shows IPV10 and IPV20. The figures 
also include the ground-based measurements from NW95 as 
horizontal lines, unless the values measured in the ground- 
based aperture fall outside the range of the vertical axes. Our 
measurements in 0.2" and 1" apertures, along with the NW95 
data, are listed in Tables[0and[0 For nearly all of the galax- 
ies there are substantial differences between the values in the 
nuclear region and the maximum STIS aperture size, as well 
as between the maximum STIS aperture size and the NW95 
measurements. These differences between the STIS apertures 
and the ground-based measurements could be due to uncertain 
resolution corrections in the narrowest lines, although galax- 
ies with narrow emission lines are not systematically more 
different than galaxies with well-resolved lines. An alternate 
explanation is that these measurements fall below those from 
NW95 because the STIS slit only subtends a fraction of the 
NLR. We will discuss this point further in Section|6] 

Figure [5] shows that the line profile widths almost always 
increase or remain approximately constant. The FWHM for 
ten galaxies increases by greater than 10%, while for an ad- 
ditional 1 1 the FWHM changes by less than 10%. Only three 
galaxies decrease by greater than 10% from the 0.2" to 1" 
aperture (NGC 1358, NGC 5427, and NGC 2110 in both 
emission lines). Two galaxies (NGC 1358 [S II], NGC 2110 
[O III]) decrease by more than 30% and are interesting cases 
because a decrease requires both significant emission at larger 
spatial scales and a velocity width that is substantially nar- 
rower at larger scales than in the nuclear region. Comparison 
of the 0.2" and 1" apertures clearly indicates that the profiles 
are substantially broader in the central region. 

IPV variations with aperture size are particularly sensitive 
to changes in the base and wings of the lines, which may 
probe the acceleration or deceleration of winds. For example, 
an increase corresponds to more high-velocity emission-line 
gas outside of the nuclear region than in the nucleus. While 
eleven of these galaxies have IPV20 variations of less than 
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FI G. 4. — IP V10 (crosses) and IPV20 (diamonds) measurements obtained from the Gaussian fits to the line profiles plotted against aperture size as in Figurel3l 
The Nelson & Whittle 1 1995) measurements for IPV10 (dotted) and IPV20 (dashed) are plotted as horizontal lines. NGC 7682 is not shown because the red wing 
of the [S II] line did not fall on the STIS CCD. 



10% between the 0.2" and 1" STIS apertures, several ex- 
hibit substantial variations. IPV20 increases by over 30% 
from 0.2" to 1" in both MKN 573 and NGC 3227, while 
it decreases by over 30% for NGC 2110 [O III] (and nearly 
this amount in [S II]). In contrast, our FW20 measurements 
for these galaxies only increase moderately with aperture 
size, if at all, which implies that the IPV width measure- 
ments are indeed more sensitive to behavior in the wings, 



as was suggested by [Whittle! l ll985al) . The substantial IPV 
decreases with radius in NGC 1358 and NGC 2110 provide 
good evidence for radiation-driven winds that decelerate at 
larger scales , as ha s already been noted for NGC 4151 by 
Nelson et al. (2000). 

We note that several IPV profiles in Figure 0] exhibit pro- 
nounced jumps in IPV value with radius between neighboring 
data values. These jumps are either due to instances where our 
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Gaussian fitting routine switched between a single-component 
fit and a multiple-component fit (or vice versa), or the pres- 
ence of knots of line-emitting gas outside of the nuclear re- 
gion. 

5.2. Asymmetries 

Figures [3] and [4] described above indicate that markedly 
different velocity components contribute to the emission-line 
profiles in the nucleus and at larger scales in approximately 
half of the sample. In addition to this information on the 
widths of these velocity components, the presence or absence 
of asymmetries can provide information on the origin of the 
line-emitting material along the line of sight, particularly in 
the presence of significant gas and dust in the NLR. 

As noted previously, the [O III] line is quite often re- 
ported to have significant, typically blue, asymmetries, 
while asymme tries are rarely observed in lower-exc itation 
lines jHeckman et alJl!981t IdeRobertis & Osterbrocklll98 
IWhittldll985bl) . Our measurements of the lWhittld ( 11985 
asymmetry parameters listed in Table ^ confirm that asym- 
metries are more common in [O III] and these values are il- 
llustrated in Figure [5] Within the 0.2" aperture we find that 
6/7 galaxies with [O III] observations have |A2o| > 0.1, while 
9/18 galaxies with [S II] observations have such asymmetries. 
While we observe asymmetries less frequently in [S II] than 
[O III], we still measure asymmetries in a substantially larger 
fraction of galaxies than typically observed in ground-based 
observations and we discuss this point further below. Two par- 
ticularly extreme cases of asymmetries are NGC 1667 [S II] 



and NGC 5194 [O III]. The asymmetries are less pronounced 
in the larger 1" aperture, where only 7/18 galaxies in the [S II] 
sample and 5/7 in the [O III] sample have such significant 
asymmetries. 

Several of these galaxies exhibit rather unusual asymme- 
tries. While most have blue asymmetries, several have sig- 
nificant red asymmetries. One peculiar case is MKN 573, 
which is more asymmetric on larger scales and this larger- 
scale asymmetry is red: A20 = -0.22 in the 1" aperture. In- 
spection of HSr images of MKN 573 Martini et alJ2003l) re- 
veals that the strong red asymmetry outside of the nucleus is 
due to the chance intersection of the STIS slit with an indi- 
vidual NLR cloud. MKN 270 [S II], NGC 1667 [S II], and 
NGC 5194 [O III] also exhibit significant red asymmetries in 
the 0.2" aperture spectra, although are less asymmetric on 
larger scales. The red-asymmetric component in these cases 
is broader than the core. 

Blue asymmetries in [O III] were suggested by 
iHeckman et al.l (119811) to be due to a combination of ra- 
dial outflow from the nuclear region and gas and dust that 
obscures the redshifted emission from the far side of the 
galaxy. The blue asymmetries confined to the nuclear region 
in NGC 4051, NGC 5347, and MKN 348 agree well with this 
interpretation. NGC 7682 and MKN 1066 havepronounced 
asymmetries on larger scales as well (see FigureQ. For these 
two galaxies 1" corresponds to several hundred parsecs and 
would require a large obscuring medium. 

We have also investigated the scale-dependence of asym- 
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TABLE 3 
IPV Width Measurements 



— Nelson & Whittle (1995) — — Aperture = 0.2" — — Aperture = 1" — — IPV20 — 

Name Line IPV20 IPV10 Line IPV20 IPV 10 IPV20 IPV 10 l"/0.2" NW95/1" 

(kms" 1 ) (kins- 1 ) (kms" 1 ) (kms" 1 ) (kms" 1 ) (kms^ 1 ) (%) (%) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



MKN 270 


[0 III] 


472 


637 


[S II] 


238 


371 


250 


393 


5.0 


89 


MKN 573 


[S III] 


423 


572 


[S II] 


132 


211 


196 


343 


48 


120 


MKN 686 


[O III] 


454 


593 


[S II] 


153 


273 


192 


333 


25 


140 


NGC 788 


[O III] 


225 


305 


[S II] 


104 


182 


117 


191 


13 


92 


NGC 1052 


[O III] 


810 


1075 


[S II] 


700 


961 


584 


829 


-17 


39 


NGC 1358 


[O III] 


355 


505 


[S II] 


339 


477 


254 


383 


-25 


40 


NGC 1667 


[O III] 


365 


525 


[S II] 


549 


714 


423 


549 


-23 


-14 


NGC 2110 


[S III] 


655 


950 


[S II] 


766 


1010 


673 


907 


-12 


-3 


NGC 2273 


[S III] 


254 


325 


[S II] 


160 


288 


189 


303 


18 


34 


NGC 3031 


[O III] 






[S II] 


441 


567 


416 


537 


-5.7 




NGC 3227 


[S III] 


775 


994 


[S II] 


394 


526 


546 


715 


39 


42 


NGC 3516 


[O III] 


490 


665 


[S II] 


202 


258 


221 


282 


9.4 


120 


NGC 3982 


[O III] 


452 


452 


[S II] 


189 


324 


198 


335 


4.8 


130 


NGC 4051 


[S III] 


293 


405 


[S II] 


224 


293 


273 


351 


22 


7.3 


NGC 4579 


[O III] 


920 


1143 


[S II] 


586 


750 


547 


704 


-6.7 


68 


NGC 5347 


[O III] 


513 


679 


[S II] 


236 


392 


233 


395 


-1.3 


120 


NGC 5427 


[O III] 


659 


956 


[S II] 


299 


554 


278 


518 


-7.0 


140 


NGC 7682 


[S III] 


266 


330 


[S II] 


774 


1102 


784 


1136 


1.3 


-66 


MKN 348 


[O III] 


543 


770 


[O III] 


669 


991 


625 


920 


-6.6 


-13 


MKN 1066 


[O III] 


640 


902 


[O III] 


445 


596 


545 


725 


22 


17 


NGC 2110 


[O III] 


655 


950 


[O III] 


1670 


2199 


1105 


1640 


-33 


-41 


NGC 4151 


[S III] 


536 


774 


[O III] 


441 


721 


533 


874 


21 


0.6 


NGC 5194 


[S III] 






[O III] 


299 


422 


232 


362 


-22 




NGC 5929 


[O III] 


412 


514 


[O III] 


254 


409 


261 


397 


2.8 


58 


NGC 7674 


[O III] 


1255 


1640 


[O III] 


1524 


1836 


1444 


1784 


-5.2 


-13 



NOTE. — Col. (1): Galaxy name. Col. (2) Line used for width measurements in NW95. Cols. (3)-(4): IPV20 and IPV10 values given by NW95. 
Col. (5): Line used for our width measurements. Cols. (6)-(7): IPV20 and IPV 10 values for our sample with a 0."2 aperture. Cols. (8)-(9): IPV20 and 
IPV 10 values for our sample with a 1" aperture. Col. (10): The percent change in IPV20 from 0."2 to 1" aperture sizes. Col. (11): The percent difference 
between our IPV20 measurement within 1" and the NW95 IPV20 value. 



metries with the spectra shown in the rightmost panels of Fig- 
ure [0 These spectra are the difference between the spectrum 
of each galaxy summed with a 1" aperture and a 0.2" aperture 
and therefore effectively isolate the emission-line component 
due to material outside of the nuclear region. The two vertical 
lines in each panel for a given galaxy correspond to the line 
peak of the emission lines in the 1" aperture. The line profiles 
in the rightmost panels relative to these vertical lines there- 
fore illustrate the extent to which emitting material outside of 
0.2" contributes to any asymmetries in the line profiles. To 
characterize the line asymmetry outside of the nuclear, 0.2" 
aperture, we calculated the fraction of the total l"-0.2" line 
flux on the blue side of the line centroid in the 1" aperture. 
A fraction greater than 0.5 corresponds to a blue asymmetry, 
while less than 0.5 corresponds to a red asymmetry. These 
values are listed in the last column of Tabled Eight galaxies 
have strong blue asymmetries outside of 0.2" (fraction > 0.6), 
two have strong red asymmetries (fraction < 0.4), and the re- 
maining 14 are either relatively symmetric (ten) or have little 
flux outside of 0.2" (five) in the STIS aperture. The signifi- 
cant asymmetry in NGC 21 10 is only seen in [O III]. 

6. DISCUSSION 

The changes in the NLR velocity field in the STIS aper- 
ture have revealed pronounced differences between the unre- 
solved, nuclear kinematics and the NLR on larger scales, yet 
these larger-scale kinematics still subtend only a fraction of 
the NLR for many galaxies because of the narrow STIS slit. 
Even the line characteristics measured from the integrated 



STIS slit may therefore differ from ground-based flux mea- 
surements. We estimated the fraction of the NLR observed 
with these observations through a comparison of the total flux 
in the STIS aperture and the valu e report ed in the ground- 
based measurements from Whittle ( 1992a). This comparison 
demonstrated that the STIS aperture inc ludes between 2 5% 
and 90% of the [O III] flux measured bv lWhittld l ll992al) for 
the seven galaxies with STIS [O III] measurements. The ve- 
locity field sampled by the STIS apertures also may depend 
on the orientation of the spectroscopic slit relative to the ma- 
jor axis of the host galaxies, as there is ev idence that rota tion 
is partly responsible for the NLR widths ( Whittldl fl 992bl) . or 
the orientation relative to the NLR and radio jet axis, which 
are known to be unrelated to the host g alaxy semimajor axis 
dUlvestad & Wilsonll984tlSchmitt et al.l2003l) . However, the 
orientation for most of these observations were not specified 
in order to reduce scheduling constraints. 

6.1. Comparison with Ground-based Measurements 

One striking characteristic of Figures[3]and|3]is how poorly 
the line profile measurements agree with the ground-based 
values on even the largest scales (see also Figure |6j. Of 
the 24 galaxies with measured FWHM in the 1" aperture, 
the FWHM of only six (NGC 1358, NGC 3031, NGC 3982, 
NGC 4051, NGC 4151, and NGC 7674) are within 10% 
of the ground-based value (see Tabled, while 12 are dis- 
crepant by greater than 30% - comparable to the observed 
scatter in the a s - cr* correlation. Surprisingly, one of these 
12 measurements (NGC 21 10 [S II]) is actually > 30% larger 



SPATIALLY RESOLVED GAS KINEMATICS IN AGN 



15 



TABLE 4 

Asymmetry Measurements 



Name 


Line 


A 1 A i"A >"1\ 

A1U (U. 2) 


A OA /A "1\ 

A20 (0. 2) 


A10 (1") 




A20 (1") 


Blue Area (1 -0. 2) 






(km s ') 


(km s ') 


(km s ') 


(km s ') 


(%) 


(1) 


& 


(3) 


(4) 


(5) 


(6) 


(7) 


MKN 270 


rs TTl 

iij 


-0.17 


-0.20 


-0.12 


-0.14 


0.43 


MKN 573 


rs TTl 


0.00 


-0.01 


-0.21 


-0.22 


0.24 


MKN 686 


rs TTl 


-0.3 


-0.03 


0.00 


0.00 




Nnr irr 


rs TTl 


^3 


11 


03 


03 


0.76 


Nnr ins? 


rs TTl 

r^ 11 1 


-0 09 


-0 09 


-0 05 


-0 06 


0.35 


. > v TV. J JJO 


rs TTl 
r^ 11 1 


01 


03 


01 


01 




NGC 1667 


[S II] 


-0.23 


-0.21 


0.03 


0.04 


0.46 


NGC2110 


[S II] 


-0.09 


-0.10 


0.00 


0.00 


0.48 


Nnr 1111 


rs TTl 


34 


0.42 


34 


39 


0.76 


Nnr ^rni 


rs tti 


00 


01 


00 


-0 01 


77 


Nnr mi 


rs TTl 


09 


12 


10 


13 


1.00 


Nnr ^si6 


rs TTl 


00 


01 


00 


01 




Nnr 398? 


rs TTl 


0.26 


0.25 


0.24 


0.22 




Nnr 4.0^1 


rs TTl 


29 


31 


25 


25 


0.82 


NGC 4579 


[S II] 


-0.06 


-0.07 


-0.08 


-0.09 


0.59 


NGC 5347 


[S II] 


0.35 


0.37 


0.32 


0.34 


0.56 


NGC 5427 


[SII] 


-0.02 


-0.02 


-0.02 


-0.02 




NGC 7682 


[S II] 










0.52 


MKN 348 


[O III] 


0.36 


0.38 


0.28 


0.29 


0.56 


MKN 1066 


[O III] 


0.41 


0.33 


0.31 


0.29 


0.46 


NGC 21 10 


[O III] 


0.15 


0.12 


0.11 


0.14 


0.78 


NGC 4151 


[O III] 


-0.02 


-0.02 


0.14 


0.12 


0.42 


NGC 5194 


[O III] 


-0.13 


-0.18 


0.03 


0.06 


0.52 


NGC 5929 


[O III] 


-0.13 


-0.17 


-0.08 


-0.06 


0.70 


NGC 7674 


[O III] 


0.46 


0.44 


0.45 


0.45 


0.75 



NOTE. — Col. (1): Galaxy name. Col. (2): Emission line used in STIS data analysis. Col. (3) - (4): A20 
and A10 asymmetry measurement with an 0.2" aperture. Col. (5) - (6): A20 and A10 measurement within a 
1 ' aperture. Col. (7): Fraction of the l"-0."2 profile area that lies blueward of the 1" line centroid. 



in the STIS aperture than in the ground-based measurement. 
This broad, nuclear component is quite obvious in Figure fj] 
For NGC 3516 the difference between the STIS and NW95 
measurements can be attributed to the slit orient ation as this 
galaxy is known to have a large, extended NLR (Pog ge et all 
1989a), yet we detect little emission in our slit outside of the 
nuclear region. 

For the seven galaxies with [O III] measurements, we exam- 
ined the difference between the FWHM in our 1 " aperture and 
the ground-based value as a function of the fraction of the flux 
within the 1" STIS aperture. The STIS flux measurement for 
N GC 7674 include s 90% of the ground-based value reported 
bv lWhittld dl992al) and our measured FWHM is within 10% 
of his value. This confirms that when most of the ground- 
based flux falls within our much narrower aperture we mea- 
sure the same kinematics. For the other six galaxies these 
data miss a larger fraction of the total flux and the line widths 
are quite different in some instances. The most notable is the 
STIS measurement of MKN 348, which includes 60% of the 
ground-based value and is 60% narrower. The STIS mea- 
surements of th e remaining fiv e galaxies include only 30% 
of the flux in the Whittle ( 1992a) aperture and the differences 
in width between the two aperture sizes range between none 
(NGC 4151) and 150% (NGC 5929). 

The disagreement between widths measured in our 1" aper- 
ture and the ground-based aperture is more striking for IPV20 
(see also Table ^ than for line widths as only two galaxies 
(NGC 4051 and NGC 4151) have STIS 1" and NW95 mea- 
surements that agree within 10%. The majority (15/24) have 
NW95 values > 30% larger than the STIS measurements and 



2/24 (NGC 7682 and NGC 2110 [O III]) are less than 30% 
of the ground-based measurements. However, we note that 
most differences in IPV20 between the STIS 1" and NW95 
are comparisons between STIS [S II] and NW95 [O III] mea- 
surements and they may therefore partly reflect more pro- 
nounced wings in the [O III] line relative to [S II]. For ex- 
ample, NGC 2110 has a 1" IPV20 = 423 km s" 1 for [S II] and 
IPV20 = 1 105 km s" 1 for [O III]. NW95 measured 655 km s" 1 
for [O III]. 

As described in Section 15.21 there are also substantial 
differences between the fraction of AGN with asymmetries 
in the low-excitation [S II] line in these data and ground- 
based measurements. This change may be due to the dilu- 
tion of nuclear asymmetries by larger-scale emission from 
host galaxy starlight or more symmetric emission from the 
NLR on larger scales. Host galaxy dilution may in particu- 
lar explain why we detect significant asymmetries in such a 
large fraction of the [S II] profiles compared to expectations 
from ground-based programs (iFilippenko & Halpernl 1 1984t 
Ide Robertis & O sterbroc kai984l:lGreene & Hol2005l) . This re- 
sult implies that the asymmetries are primarily due to nuclear 
line-emitting gas, rather than material more evenly distributed 
throughout the NLR on larger scales. 

This comparison of width and IPV values suggest that mea- 
surements within the STIS aperture are generally smaller than 
the ground-based values. A quantitative comparison demon- 
strates that this is the case, relative to both ground-based 
NLR width measurements and er* (see Tableland Figure^. 
Comparison of the STIS FWHM measurements in 0.2", 1", 
lOOpc, and 200pc apertures (divided by 2.354 to approximate 
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FIG. 5. — A20 measurements for the [S II] and [O III] emission lines in our sample at aperture sizes of 0.2" (crosses) and 1" (diamonds). Positive A20 values 
indicate blue asymmetries, negative A20 values indicate red asymmetries, and near-zero values indicate symmetric profiles. The average A20 values for [S II] 
and [O III] at each aperture are shown. 



a Gaussian <j g ) to er* show that the STIS measurements sys- 
tematically underestimate cr* by 10 - 20%. The STIS mea- 
surements similarly underestimate the ground-based NW95 
FWHM measurements. Figure [6] does not show any signifi- 
cant evidence for differences between [S II] and [O III]. The 
scatter between STIS and cr* measurements is slightly worse 
(30 - 40%) than the scatter between STIS and NW95 FWHM 
measurements (20 - 40%). In contrast, FigureQindicates that 
STIS measurements at 0.2" and 1" and lOOpc and 200pc are 
significantly more similar in value and have substantially less 
scatter (10 - 20%, see also Table^. The smaller measured 
widths in the STIS slit is likely due to the collisional nature 
of gas, which will therefore tend to reside at least partially in 
a disk. A narrow slit at a random orientation will then sample 
a smaller fraction of the full radial velocity field than a larger 
aperture slit and therefore measure a smaller width. 

The comparable number of galaxies with red and blue 
asymmetries on larger scales, and more importantly that these 
asymmetries are in general substantially weaker than those in 
the 0.2" aperture, suggest that asymmetries measured in large 
(including ground-based) apertures originate on small scales, 
or less than ~ lOOpc based on the spatial resolution of these 
observations. 

6.2. Correlations with Global Properties 

In addition to the quality of the correlations between 
line widths measured in various apertures and er*, we also 
searched for systematic trends between the residuals and the 
properties of the galaxies. We specifically investigated the 
residuals between the measurements in the 1" aperture and cr» 



as a function of distance, Hubble type, fraction of the galaxy 
in a 1" aperture, and radio power. Significant residuals as a 
function of galaxy distance, Hubble type, and fraction of the 
angular size of the galaxy subtended by the spectroscopic slit 
would indicate very useful information about the origin of the 
scatter in the o^-a* relation and these quantities could poten- 
tially be used to empirically determine corrections to reduce 
scatter, as well as provide insight into its origin. For example, 
the NLR measured exactly within the host bulge's effective 
radius might prove to be the best tracer of ov Some evidence 
for such a relation might be revealed in the residuals plotted 
against these three parameters, however as shown in Figure[8] 
there is no evidence that the residuals are correlated with any 
of thes e parameter s. 

Whittle ( 1992c) found evidence that Seyfert galaxies with 
linear radio morphology and high luminosity tend to have 
broader lines. We do not see significant correlation between 
residuals in a g - er» in our data, although only three galax- 
ies in our sample are above his radio luminosity threshold 
of Lhi5 > 10 W Hz" 1 . Other parameters mentioned in 
the Introduction that correlate with sy stematically b roader 
emission-line widths include L/LehiI ifpreene & Hoi 120051) 
and disturbed morphologies (Whittle 1992b), however this 
sample does not contain a sufficient number of objects with ei- 
ther high accretion-rates or significantly-disturbed morpholo- 
gies to investigate these quantities. Boroson (2005) also re- 
cently noted that strongly blueshifted [O III] lines (relative to 
the host galaxy) were systematically broader than similar ob- 
jects that were not blueshifted. These objects tended to be 
those with large L/LEdd- None of the seven galaxies in our 
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FIG. 6.— Our measurements of FWHM/2.354 plotted against NW95 cr, (left panels) and NW95 FWHM/2.354 (right panels). The top panels (a & b) show our 
[S II] (diamond) and [O III] (star) FWHM measurements in an 0.2" aperture, the middle two panels show our measurements in a 1" aperture, and the bottom two 
panels show our measurements in a fixed physical aperture size of lOOpc. The dashed line in each frame represents a 1:1 correspondence between the quantities 
plotted. 
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TABLE 5 
Statistical Results 





Aperture Size 


[S II] 


[O III] 


Total 






(kms- 1 ) 


(kms -1 ) 


(kms- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


((FWHM/2.354)/ct„) 


0."2 


0.81 ±0.36 


0.76 ±0.11 


0.80 ± 0.34 




1" 


0.82 ±0.31 


0.80 ± 0.33 


0.81 ±0.31 




lOOpc 


0.83 ±0.31 


0.81 ± 0.27 


0.83 ± 0.30 




200pc 


0.85 ± 0.30 


0.83 ± 0.30 


0.84 ± 0.29 


(FWHM/NW95 FWHM) 


0."2 


0.84 ± 0.33 


0.70 ± 0.34 


0.80 ± 0.33 




1" 


0.84 ± 0.20 


0.70 ± 0.22 


0.80 ± 0.21 




lOOpc 


0.86 ± 0.27 


0.72 ± 0.35 


0.82 ± 0.29 




200pc 


0.88 ± 0.24 


0.70 ± 0.21 


0.83 ± 0.24 


(FWHM/FWHM) 


l"/0."2 


1.06 ± 0.26 


1.07 ± 0.23 


1.06 ±0.25 




200pc / lOOpc 


1.03 ±0.10 


1.04 ± 0.20 


1.03 ±0.13 


((NW95 FWHM/2.354)/o-» ) 






1.15 ±0.33 





NOTE. — Col. (1): Pair of values used for statistical comparison. Col. (2): Aperture size 
used in our measurement of FWHM/2.354. Col. (3): Mean and standard deviation of the 
sources for which we measure the [S II] line widths. Col. (4): Mean and standard deviation 
of the sources for which we measure the [O III] line widths. Col. (5): Mean and standard 
deviation of all sources. 



300 



200 

in 
E 



100 



O [Sll] FWHM/2.35 
* [ONI] FWHM/2.35 



O 



' o 



V' 

o / 

"f&O o 



'O 



o 



(a) 



100 200 300 

FWHM/2.35 (km s"') in 0."2 



o 
o 



(0 

£ 

fO 
(N 
\ 

X 




100 200 300 

FWHM/2.35 (km s"') in 100pc 



FIG. 7. — Comparison of our [S II] (diamonds) and [O III] (stars) FWHM measurements in different apertures, (a) Width measured in a 1" aperture plotted 
against the width measured in an 0.2" aperture, (b) Width measured in a 200pc aperture plotted against the width measured in a lOOpc aperture. The dashed line 
in each frame represents a 1:1 correspondence between the quantities plotted. 
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FIG. 8. — Comparison of residuals between gas and stellar kinematics and various global properties. The four panels show: (a) Distance; (b) Hubble T type; 
(c) Fraction of the galaxy size (2" divided by the angular size provided in Table IT]; and (d) radio power. The points to the left of the dashed line do not have a 
measured radio flux and are only shown to indicate the ordinate values. 



[O III] sample have substantial blueshifts. 

7. SUMMARY 

We have conducted a detailed analysis of spatially-resolved 
[O III] and [S II] NLR emission from 24 well-studied, nearby 
AGN. These observations have detected considerable emis- 
sion outside of the unresolved nucleus (0.2" or 10 - lOOpc) 
and this emission often contributes significantly to the mea- 
sured line profiles at larger scales. We have characterized the 
spatial dependence of this emission with a range of width, 
area, and asymmetry measurements and shown there are not 
only substantial changes from 0.2" to 1" (typical of ground- 
based observations), but also from a 1" long STIS aperture 
(with either a 0.1" or 0.2" width) and the approximately 
1" x 1" or larger aperture used for ground-based measure- 
ments of these galaxies. In particular, the large variations in 
line profiles as a function of aperture size demonstrates that 
the profile of the NLR lines are set by the kinematics of the 
gas in the NLR itself, and not just radiatively-driven winds 
from the nucleus. 

The spatial scale(s) responsible for the NLR emission lines 
has received substantial recent attention because the widths 
of the NLR lines may be used to estimate the stellar veloc- 
ity dispersion 0%, which in turn can be used as a proxy for 
the black hole mass M # . Our analysis shows that even at the 
largest spatial scales observed with STIS, the line widths are 
systematically smaller by 10 - 20% than ground-based width 
measurements (as well as cr*). As we estimate that the STIS 
slit may include as little as 25% of the NLR flux, emission 



on larger scales or outside of the narrow STIS slit width must 
broaden the line profiles. 

While the line profile measurements from the STIS data are 
systematically less than the ground-based NLR and cr» values, 
the scatter is comparable. The substantial scatter in the <j g — cr* 
relation therefore appears to be due to a sufficiently complex 
set of parameters that it is effectively stochastic, at least to 
the extent that the scatter can not be reduced to the magnitude 
of the scatter in the M, - a* relation. In addition to rotation, 
Eddington ratio, and compact radio jets, the list of parame- 
ters that increase the scatter should include the dumpiness 
of the NLR, the orientation of the NLR with respect to the 
host galaxy semimajor axis, and the amount and distribution 
of dust in the NLR. If the NLR were as isotropic as the stellar 
distribution, the scatter would be comparable to the M m - cr* 
relation, yet because the gas is collisional, and also clumpy, it 
can not be as good a tracer of the bulge potential. 

This interpretation is supported by the substantial and di- 
verse types of width variations observed in these spatially- 
resolved measurements. Specifically, the most common types 
of line width variation are < 30% changes between the 0.2", 
1", and NW95 apertures (nine galaxies) and relatively little 
width increase in the STIS aperture, but a substantial increase 
(> 30%) between the 1" STIS and NW95 measurements (also 
nine galaxies). The remaining types of variations observed 
are greater than 30% width changes on the STIS scales only 
(four) and on all scales (two). Two of the four galaxies with 
only substantial profiles changes on the STIS scale have sub- 
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stantially broader emission in the 0.2" than in the 1" aper- 
ture. These large nuclear velocities may be due to radiatively- 
driven winds that decelerate or simply do not extend to larger 
scales. However, neither of these nuclear line profiles exhibit 
substantial asymmetries, which suggests that if this emission 
is due to a uniform outflow, there is not substantial dust in the 
nuclear region. 

Both red and blue asymmetries are observed in the [S II] 
and [O III] line profiles. These observations are unusual in two 
respects. First, asymmetries in [S II] are reported far less fre- 
quently than in [O III], yet we observe [S II] asymmetries with 
nearly comparable frequency. The low frequency of observed 
[S II] asymmetries in ground-based observations have been 
ascribed to the fact that [S II] is a lower-ionization line. Our 
observations suggest that asymmetries are commonly present 
in [S II], but the weaker [S II] emission is more easily di- 
luted by host galaxy starlight and not as readily observed 
in the larger ground-based apertures. Secondly, we observe 
several examples of galaxies with red asymmetries in either 
[S II] or [O III]. For the cases in which the red asymmetries 
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are largely confined to the unresolved nuclear spectrum, these 
observations do not agree well with outflow models that pro- 
duce blue asymmetries through invocation of obscuration of 
the far side of the galaxy. These red asymmetries may be due 
to the patchy nature of the ISM on very small scales, which 
may produce uneven illumination of NLR clouds and a patchy 
line-of-sight velocity distribution for the NLR clouds. Such 
variations within the NLR may be responsible for the bulk of 
the scatter in the a„ - correlation. 
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